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FOREWORD 


The  transportation  fuel  industry  is  experiencing  a  time  of  change. 
Depletion  of  fossil  petroleum  reserves  and  increased  costs  of 
petroleum-derived  fuels  are  forcing  the  industry  toward  alternative  fuels. 
There  is  also  an  urgent  need  to  develop  cleaner-burning  fuels  to  relieve  the 
pressures  of  air  pollution  on  our  environment. 

The  best  current  option  for  a  future  transportation  fuel  base  appears  to 
be  methanol,  a  simple  alcohol  that  can  be  abundantly  produced  from  natural 
gas,  coal  or  biomass.  When  synthesized  from  coal  or  biomass,  the  result  is  i 
mixture  of  methanol  with  small  quantities  of  higher  alcohols ,  such  as 
ethanol,  butanol  and  propanol.  This  methanol  plus  higher  alcohols  mixture 
(MPHA),  in  combination  with  gasoline  and  additional  ethanol,  is  likely  to 
become  the  transportation  fuel  of  the  future.   With  some  modifications, 
existing  automobile  engines  can  use  methanol  more  efficiently  than  gasoline 
and  simultaneously  release  fewer  emissions. 

This  project  developed  base-line  data  on  the  physical  and  combustion 
characteristics  of  MPHA  in  mixtures  with  gasoline.   The  data  should  prove 
invaluable  to  automotive  engineers  for  designing  more-efficient  and  cleaner 
operating  engines  capable  of  using  wide  ranges  of  fuel  mixtures.   ENR  is  . 
pleased  to  have  sponsored  this  research  and  hopes  that  the  information 
provides  a  useful  base  for  redefining  our  transportation  energy  resources. 


Don  Etchison,  Director 
Illinois  Department  of  Energy 
and  Natural  Resources 
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CHAPTER  I 


INTRODUCTION 


Methanol  received  wide  attention  as  a  fuel  for  the  spark  ignition 
engine  in  pure  form  or  as  a  blend  with  gasoline.  Methanol  today  can 
be  produced  from  the  nation's  abundant  resources,  coal,  waste 
products  from  farms  and  municipalities,  as  possible  additional 
f eedstock[l-10]*.  It  is  probable  that  the  synthetically  produced 
methanol  fuel  of  the  future  will,  in  fact,  be  a  mixture  of  methanol 
and  higher  alcohols,  rather  than  pure  methanol.  In  the  low  pressure 
methanol-product ion  process,  the  purity  of  methanol  is  achieved  by 
conventional  distillation  methods.  The  other  components  are  higher 
alcohols  which  can  be  extracted  at  lower  levels  of  the  distillation 
column  [12].  The  use  of  methanol  plus  higher  alcohols  instead  of 
pure  methanol  reduces  the  production  cost  of  the  fuel,  in  general. 


The  composition  of  alcohol  fuels  derived  from  coal  and  other 
sources  depends  upon  the  manufacturing  process,  process  operating 
conditions,  and  properties  of  the  raw  material  [11].  The  alcohol 
blend  which  is  representative  of  the  product  which  could  be  obtained 
from  a  methanol  synthesis  process  using  a  non-selective  catalyst 
contains  by  volume,  75%  methanol,  5%  ethanol,  7.5%  1-propanol,  and 
12.5%  2-methyl-l-propanol  [13].  This  composition  of  alcohols  was 
used  as  the  Methanol  Plus  Higher  Alcohols  fuel.  Because  of  different 
comocsitions  cf   indolene-methanol  blends  and   indciene-MPHA  blends, 


*  Numbers  in  brackets  refer  to  references  in  Chapter  VII. 


the  physical  properties  and  performance  parameters  need  to  be 
examined  at  various  operating  conditions.  The  fuel  specifications, 
are  listed  in  Table  1.1.  One  of  the  most  important  properties  of 
methanol  and  MPHA  is  the  high  octane  number,  which  allows  its  use  at 
higher  compression  ratios,  resulting  in  an  increase  in  the  thermal 
efficiency  of  the  spark  ignition  engine.  The  objective  of  this  work 
is  to  investigate  the  physical  and  chemical  properties  of  MPHA- 
indolene  blends,  and  the  combustion  characteristics. 

The  experiments  were  performed  at  various  methanol  and  MPHA 
concentrations  (from  0  to  100%)  in  the  blend.  The  compression  ratio 
was  varied  from  5:1  to  Knock  Limited  Compression  Ratio  (KLCR).  The 
fuel  properties,  performance  and  combustion  characteristics  were 
also  evaluated. 

In  Chapter  II,  the  experimental  facility  has  been  described.  In 
Chapters  III,  IV  and  V,  the  data  analysis,  theoretical  analysis,  and 
results  are  presented,  along  with  a  discussion.  Chapter  VI  contains 
the  conclusions  and  recommendations. 


TABLE    1.1 


Puel    In»p«ction   Data 


Property 


XadolMM    Itothaaol    Bthoaol    ■  TwuprnMOl         Imbbatmnol 


For aula 


CH3OH 


C2H3OH 


C3H7OH    CH3CH(CH3)CH20H 


Specific  gravity 
60»F/60*F 


0.75 


0.793 


0.8111 


0.8045 


0.8030 


Boiling  t«ap. 

Or 


80-430 


147 


171.1 


207 


226.2 


Flash  Point 
0.. 


-45 


54 


68 


81 


95 


Autoignition  teap.     495 


Vapor  pressure 
at  68 *F.  HMHG 
Flaaability  limits 


97.68 


48.1 


14.9 


lovar,  % 


7.6 


1.3 


36 


6.0 


3.3 


13.7 


2.1 


10.6 


1.7 


Vapor  density 
(air  -  1) 


1.1 


1.3 


2.1 


2.5 


CHAPTER  II 


EXPERIMENTAL   SETUP 


Engine  and  Instrumentation 

A  standard  single  cylinder  CFR  spark  ignition  engine  with 
variable  compression  ratio,  manufactured  by  Waukesha  Motor  Company 
was  used  in  these  experiments.  The  combustion  chamber,  with  a  flat 
cylinder  head,  is  shown  in  Fig. 2.1.  One  minor  modification  was  made 
to  facilitate  measurement  of  the  flame  speed.  A  spark  plug  with 
extended  electrodes  was  placed  in  the  magnetic  pick-up  plug 
position,  and  an  ionization  probe  was  placed  in  at  the  spark  plug 
position.  This  was  done  in  order  to  allow  deep  probing  into  the 
combustion  chamber,  without  interfering  with  the  piston  motion.  This 
engine  was  coupled  to  a  General  Electric  CD-77  dynamometer,  which 
can  absorb  30  horsepower  and  deliver  25  horsepower.  Fig. 2. 2  shows 
the  overall  layout  of  the  experimental  setup. 

A  Bendix  breakerless  ignition  system  is  used  in  this  engine.  This 
system  consists  of  an  ignition  coil,  a  control  unit,  a  spark  plug 
cable,  and  a  spark  plug.  An  AC  spark  plug  with  extended  electrodes 
was  cut  to  allign  it  with  the  center  line  of  the  travelling  probe, 
as  in  Fig. 2  .  3 


The  intake  air  system  consists  of  an  air-filter,  a  laminar  flow 
element  (Meriam  50MW20-2),  and  an  air  heater.  The  absolute  intake 
air  pressure  is  measured  at  the  intake  of  the  laminar  flow  element, 
with  an  inlet  absolute  pressure  gage  (Meriam  310EF10).   The  air  flow 
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Fig.    2.1      Top   and    Side   Views    of    C.F.R.    Combustior.   Chamber 


Instrumentation 
Systems 


Ionization 


pyoccny** 


transducer 


Trigger 


10  C.A. 
1  C.A 


Fly 
wheel 


Cylinder 
head 


Engine 


Dynamometer 


Magnetic 
pickup 


load  pickup 


Enclosure 


Speed 
variator 


Air-fuel  mixture 


Torque 
variator 


Intake 
manifold 


K) 


Compression 

ratio 

indicator 


Dry  type 
transformei 


RPM 
indicator 


LOAD 
indicator 


Swich 


Power 
source 


Fig. 2*2  Overall  engine  and  instrumentation  systems 
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Fig.  2.3  Spark  plug  and  its  alignment  with  th^  probe 


rate  is  computed  from  the  pressure  differential  across  the  laminar 
flow  element,  with  a  differential  pressure  indicator  (Meriam 
40HE35).  The  inlet  air  temperature  is  kept  constant  at  120  F,  by  an 
electric  heater.  Fig. 2. 4  shows  the  layout  of  the  fuel-air  intake 
system. 

The  engine  speed  was  measures  by  a  shaft  encoder  (BEI- 
H250-50-A-7460R-EM16) ,  and  a  frequency  counter.  The  engine  load  was 
measured  with  a  load  cell  (Interface  SSMAJ250)  and  a  digital 
transducer-indicator  (Newport  Q2000SGR3).  The  load  cell  is  used  to 
measure  the  engine  torque  and  to  find  the  minimum  spark  advance  for 
best  torque  (MBT). 

A  high  frequency  response  AVL  pressure  transducer  type  KQP300CVK 
was  used  to  measure  the  cylinder  pressure.  This  transducer  is  water 
cooled,  and  is  vibration  and  temperature-  compensated.  It  is  used 
in  conjunction  with  a  Kistler  504E  Charge  Amplifier.  Prior  to 
measuring  cylinder  pressure,  the  transducer  was  coated  with 
silicone,  and  calibrated  along  with  its  mating  amplifier. 

The  fuel  system  consists  of  a  fuel  reservoir,  a  fuel  puT.p,  a  fuel 
container,  a  weighing  scale,  an  electronic  timer,  a  fuel  fuel 
filter,  and  an  injection  nozzle.  Clear  indolene,  methanol,  ethanol, 
propanol,  and  butanol  (industrial  grade)  vere  used  in  this  study. 
The  rate  of  fuel  consumption  is  measured  by  a  weighing  scale  and 
timer. 

Flame  propagation  time  is  measured  with  the  ionization  probe.  It 
consists  of   a  water   cooled  probe,    and   an  electronic   amplifier- 
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10 


counter.  Fig. 2. 5  shows  the  layout  of  the  ionization  probe.  This 
probe  was  designed,  and  manufactured  here,  at  U.I.C.  The  probe  can 
traverse  along  its  center  line,  inside  the  combustion  chamber,  and 
can  be  held  at  any  desired  distance  from  the  spark  plug.  The 
electrodes  of  the  ionization  probe  consist  of  two  thermocouple 
wires-  a  platinum  one  and  a  platinum-13%rhodium  wire.  They  are 
separated  by  MgO  insulation,  and  enclosed  in  an  inconel  sheath.  The 
diameter  of  the  wires  is  0.063  in.  The  outer  ends  of  the  electrodes 
are  connected  to  the  electronic  amplifier-counter,  and  to  an 
oscilloscope,  the  signal  from  which  is  input  to  the  data  acquisition 
system. 

Data  acquisition  and  recording  system 

The  data  acquisition  and  recording  system  consists  of  1  CAD  at 
2000  rpm  acquisition  rate,  2  A  to  D  channels,  1  interval  timer,  1 
TDC  detector,  1  parallel  data  latch,  an  IBM  XT  PC,  and  an  FX85 
printer.  Fig. 2. 6  shows  the  layout  of  the  data  acquisition  system. 
The  following  data  were  recorded  and  monitored  on  the  oscilloscope: 
1.  Instantaneous  cylinder  pressure  signal.  2.  Spark  signal  which 
indicates  the  spark  timing.  3.  Ionization  probe  signal.  4.  Shaft 
encoder  signal. 


In  addition,  mass  flow  rate  of  air  and  fuel,  torque,  tempe  atures 
of  oil,  coolant,  and  exhaust  gas  were  also  recorded  to  obtain  the 
performace  parameters. 
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Equipment  Preparation 

The  engine  was  given  a  complete  overhaul,  including  a  cleaning  of 
the  piston,  cylinder  head,  intake  and  exhaust  valves,  intake  and 
exhaust  manifolds,  cylinder  pressure  transducer,  injection  pump  and 
cooling  system.  The  pressure  transducer  was  calibrated  and  mounted 
flush  with  the  cylinder  head.  The  load  cell-indicator  system  was 
also  calibrated.  The  intake  and  exhaust  valves  were  carefully  lapped 
to  ensure  proper  seating  before  installation  in  the  engine. 
Calculations  were  made  to  determine  the  stoichiometric  air  to  fuel 
ratios  for  all  the  fuel  blends,  as  described  in  the  appendix.  The 
data  recording  sheets  were  developed  as  seen  in  the  sample 
calculation. 

The  experiments  were  divided  into  two  parts. 

The  first  part  was  to  obtain  the  Minimum  spark  advance  for  Best 
Torque  (MBT) ,  and  the  Knock  Limited  Compression  ratio  (KLCR)  with 
indolene-methanol  and  indolene-MPHA  blends.  The  second  part  was  to 
obtain  detailed  combustion  and  performance  characteristics  of  all 
fuel  blends.  The  engine  operating  conditions  and  test  matrix  are 
given  in  Table  2.1-2.3. 

Experimental  Procedure 


The  following  experimental  procedure  was  used: 

1.  Prior   to   running   experiments,    the   engine   and  other 
equipments  were  prepared  as  stated  in  the  previous  section. 

2.  The  checklist   shown  in  the  appendix  was  followed  to  motor 
and  fire  the  engine. 
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TABLE  2.1 


TEST  MATRIX  FOR  INDOLENE -METHANOL  BLENDS 


CR 


100%I 


70%I+ 
30%M 


50%I+ 
50%M 


30%I+ 
70%M 


100%M 


5.0 


101 


104 


108 


112 


116 


6.0 


102 


105 


109 


113 


117 


7.0 


103 


106 


110 


114 


118 


8.0 


107 


119 


9.0 


120 


10.0 


111 


115 


121 


'  •  ■<>■ 
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TABLE  2.2 


TEST  MATRIX  FOR  INDOLENE-MPHA  BLENDS 


CR 


100%I 


70%I  + 
30%MPHA 


50%I+ 
50%MPHA 


30%I  + 
70%MPHA 


100%MPHA 


5.0 


101 


122 


126 


130 


134 


6.0 


7.0 


102 

103 


123 


124 


127 


128 


131 


132 


135 


136 


8.0 


9.0 


125 


10.0 


129 


133 
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Table    2.3 
GENERAL  OPERATING  CONDITIONS 


Spark  plug   gap   (in.) 

Spark   timing 

Valve   clearance    (in.) 

Oil    (crankcase) 

Intake   air    temperature    ( 'F) 

Equivalence    ratio 

Speed    (RPM) 

Compression   ratio 

Fuel 


Injection   pressure    (psi) 
(manifold 


0.045 

MBT 

0.008    +   0.001 
(hot) 

Synthetic  oil 
design    for  methanol 

120    +   2 

1    +   0.02 

1000    +   5 

5    to   KLCR 

indole ne  ,    indolene- 
alcohols     blends. 

1200 


Intake    air   pressure 


WOT 
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3.  The  compression  ratio  was  adjusted  and  varied  to  obtain  the 
KLCR. 

4.  The  engine  speed  was  set  at  a  constant  value  of  about  1000 
rpm. 

5.  The  equivalence  ratio  was  set  to  about  1.0  by  using  the 
calculation  sheet  shown  in  the  appendix. 

6.  The  engine  was  run  for  a  period  of  an  hour,  at  the  operating 
conditions  given  in  Table  2.3  to  ensure  steady  state 
operation  before  recording  the  data. 

7.  In  the  first  part,  the  spark  advance  was  changed  from  20 
deg.BTDC  to  50  deg.BTDC  in  increments  of  one  degree  and  the 
load  cell  reading  was  recorded.  The  spark  advance  at  which 
the  load  cell  indicated  maximum  torque  is  called  the  MBT 
spark  advance. 

8.  The  compression  ratio  was  varied  from  5:1  to  KLCR  in 
increments  of  1. 


Performance  Parameters 

The  parameters  which  have  been  studied  in  this  work  are;  1.  brake 
power,  2.  brake  specific  fuel  consumption  and  3.  brake  thermal 
efficiency.  The  tests  started  with  pure  indolene  .  Methanol  was  then 
;idded  by  30%,  50%  and  70%.  Tables  2.1  and  2.2  show  the  experimental 
test  matrix.  The  compression  ratio  was  varied  from  5:1  to  KLCR, 
i.e.,  when  engine  knock  became  audible.  The  knock  was  observed  from 
the  pressure  trace. 

Burning  veloc  i ty 

The  parameters  which  have  been  studied  in  this  part  of  the  study 
are:  1. Flame  speed,  and  2. Burning  velocity.  The  volume  fraction  of 
methanol  and  MPHA  were  varied  by  0.1,  0.2,  0.3,  and  0.4.  The 
compression  ratio  was  kept  constant  at  5:1.  The  flame  propagation 
time  was  measured  at  eight  locations  during  flame  propagation. 
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CHAPTER  III:  DATA  ANALYSIS  AND  SAMPLE  CALCULATION 


In   this   chapter,    the  data  analysis  and  a   sample  of   the 

calculation  will   be  presented.   The  data  analysis  consists  of  two 

parts:    1.   Flame  travel  time  analysis.    2.   Pressure-Time  history 
analysis. 

Flame  Travel  Time  Analysis 

The  flame  travelling  time  is  taken  as  the  time  interval  between 
the  spark  timing  and  the  arrival  of  the  flame  at  the  tip  of  the 
ionization  probe.  A  high  speed  time  counter  is  used  to  measure  the 
flame  travelling  time  in  milliseconds.  The  spark  and  ionization 
signals  are  also  monitored  on  the  oscilloscope.  For  each  run,  512 
data  points  are  temporarily  stored  in  the  microcomputer  memory  and 
later  transferred  to  files. 

After  determining  the  512  data  points  for  the  flame  propagation 
time  for  each  position  of  the  ionization  probe,  an  arithmetic 
average  is  performed  on  the  data  by  the  program  TIMEAVG.  This 
program  also  checks  and  compares  the  mean  with  all  the  data  points, 
excludes  bad  data  points,  and  recalculates  the  average  flame  travel 
time  (Bad  data  "Jints  were  defined  as  above  1.5  times  or  below  0.5 
times  the  mean) . 

This  process  was  repeated  for  all  the  different  probe  positions 
between  the  spark  plug  and  the  chamber  wall.   The  time  average  data 


were  then  curve  fitted  using  a  fifth  degree  polynomial  to  give  the 
distance  vs.  time  relationship  for  the  flame  travel  process.  In 
order  to  improve  the  accuracy  of  the  curve  fitting  routine, 
additional  points  were  inserted  near  the  end  points. 

The   apparent  flame   speed   is  computed   from   the  slope  of   the 

distance-time  curve   by  a   computer  program  APPBVEL.  This   program 

calculates  the  apparent  flame  speed  at  every  degree  crank  angle  by 
using  the  polynomial  coefficients  as  input. 


Pressure-Time  History  Analysis 

The  cycle-to-cycle  variation  in  pressure  traces  is  a  common 
problem  in  spark  ignition  engines.  The  combustion  characteristics 
calculated  from  the  model  are  very  much  dependant  upon  the  pressure, 
and  its  derivatives.  An  attempt  has  been  made  in  the  literature  to 
use  the  average  pressure  of  many  cycles.  Pattersons [18 ]  has 
recommended  that  2000  pressure  cycles  be  used  for  a  representative 
average  cycle,  and  Peter[19]  has  recommended  300  cycles.  Recently 
Mattavi  et  al  [20]  used  96  consecutive  cycles  to  obtain  the  average 
pressure  data.  Following  was  the  procedure  used  in  this  study  to 
obtain  the  average  pressure  data.  The  high  speed  data  acquisition 
system  was  used  to  collect  80  consecutive  cycles  for  each  run.  The 
pressure  was  recorded  at  every  1/2  deg.  crank  angle.  The  arithmetic 
average  of  these  cycles  was  taken  to  obtain  one  average  pressure 
cycle. 
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Sample  Calculation 


'M\ 


Flame  Speed 

The  average   time  of   512  points   for  eight  different   probe 
positions  are  shown  below: 


Flame  Travel  Distance 
Time 

0.85000E-02  8.49 

0.80000E-02  8.50 

0.75000E-02  8.25 

0.72186E-02  8.09 

0.62738E-02  7.09 

0.56755E-02  6.09 

0.49691E-02  5.09 

0.42186E-02  4.09 

0.36305E-02  3.09 

0.28760E-02  2.09 

0.22500E-02  1.50 

0.17868E-02  1.09 

0.10377E-02  0.69 

0.07500E-02  0.50 

0.05000E-02  0.25 

0.02500E-02  0.15 

O.OOOOOE-02  0.00 


Improved  Flame  Travel 
Time 


0.72186E-02 
0.64068E-02 
0.56755E-02 
0.49691E-02 
0.43507E-02 
0.36305E-02 
0.28760E-02 

0.17934E-02 
0.10979E-02 


Fig.   3.1  is  a  plot  of  the  distance  vs.   flame  travel  time.   The 

J 


circular  symbols  are  the  inserted  data  points  to  increase  the 
accuracy  of  the  polynomial  curve  fitting.  After  curve  fitting,  the 
coefficients  are  used  to  calculate  the  flame  speed  at  every  crank 
angle.  The  appendix  contains  plots  of  flame  speed  vs.  crank  angle 
calculated  from  APPBVEL  program. 


^ 


Mass  Fraction  Burne<' 


As  shown  in  Fig.  3.2,  ignition  delay  and  combustion  interval  are 
calculated  from  the  mass  fraction  vs.  degree  crank  angle  curve. 
These  data  are  obtained  from  the  computer  simulation  program  SPMIX. 
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DISTRNCE  VS  FLAME  TRRVEL  TIME  AT  CR-5.0t l.INDOLENE 


0. '5  0.30  0.  A5  o.ao 

TI)1E^lD*^+2  (SECOND) 


0.75 


0.  90 


5   i^StS'3 


QRTR 


Fig.    3.1        Distance-Time   curve 
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MASS  FRACTION  BURNED  VS  DEGREE  Cfl 


-30.  00 


-10.00     0.00 
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20.00 


MfiSS  FRfiCT,   BURNED 


Fiq.    3.2      Mass    fraction   burned  vs    **CA 


Determinat  ion  of  Combustion  Characteristics 

The  fundamental  quantities  from  which  combustion  can  be 
characterized  are,  actual  burning  velocity,  ignition  delay,  and 
combustion  interval.  The  following  is  an  outline  of  the  procedure 
used  in  calculating  these  fundamental  quantities. 

Ignition  Delay 

The  period  of  time  from  spark  to  the  point  when  heat  release  is 
detected,  is  called  ignition  delay.  This  period  can  be  calculated 
from  the  p-t  diagram  where  the  pressure  increases  drastically* 
Alternatively,  this  period  can  be  obtained  from  the  mass  fraction 
burned  vs.  degree  crank  angle  curve.  Ignition  delay  can  be  defined 
as  the  interval  between  the  spark  and  the  time  when  5%  of  the  mass 
is  burned.  This  definition  is  used  in  this  study  to  calculate  the 
ignition  delay  period.  From  Fig.  3.2,  it  can  be  seen  clearly  that 
after  5%  of  the  mass  is  burned,  the  rate  of  burning  rises  sharply. 

Combustion  Interval 


After  the  ignition  delay  period,  the  fully  developed  turbulent 
flame  propagates  across  the  combustion  chamber.  Combustion  interval 
is  the  period  of  time  from  the  end  of  the  ignition  delay  period  to 
the  end  of  combustion.  Keeping  in  mind  the  accuracy  of  the  model, 
the  combustion  interval  is  defined  as  the  time  between  5%  to  90% 
mass  burned.  During  this  interval,  the  flame  propagates  with  a 
finite  thickness  through  the  unburned  mixture.  Fig.  3.2  shows  the 
definition  of  combustion  interval  from  the  mass  fraction  burned  vs. 
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degree  crank  angle  curve. 

Burning  Velocity 

The  burning  velocity  can  be  calculated  from  the  flame  speed  and 
the  pressure-time  data,  the  flame  front  area  at  every  crank  angle  is 
needed  for  these  calculations.  Flame  front  area  was  calculated  from 
the  spherical  flame  propagation  model.  Fig.  3.4  is  a  plot  of  the 
flame  front  area  vs.  flame  radius. 

Fig.  3.5  shows  the  curves  of  flame  speed,  expansion  velocity,  and 
actual  burning  velocity  as  functions  of  crank  angle  degree,  for 
indolene. 


MASS  AND  VOL.  FRflCT,  BURNED  VS  DEGREE  Cfl 


lV-'/  ■ 

^;:^^■ 


-30.  00 


•20.00     -10.00     0.00 

DEGREE  CRANK  ANGLE 


fO.  00 


20.00 


^•■■.'■■: 


t 


MRSS  FRflCT.   BURNEJ 

VOL.  FPK^r.  SSSS 


Fig.  3.3    Plot  of  mass  and  volume  fraction  burned  vs  degree 
crank  angle  at  base  conditions. 
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FLflME  RffilUS  VS  SPHERICPL  FLflME  FRONT  fiRER  AT  CR-5 
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Fig.  3.4  Flame  front  area  vs  flame  radius 
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VELOCITY  VS  DEGREE  Cfl  FOR  INDOLENE,  CR=5 
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Fig,  3.5   Burning  Velocity  vs  'CA 
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CHAPTER  IV 


THEORETICAL  ANALYSIS 


^''^      Models  of  Combust  inn 

MoaeUn,  of  t.e  cc.ustion  phenomena  .ns.de  the  spar, 
xgnition  engine  is  relatively  si^pie  if   compared  „ith  .ode- 
ling  of  diesel  combustion. 

in  spar,  ignition  engine,  the  fixture  of  air  and  f.el 
"  -— P-i-a  ana  can  .e  assumed  as  an  homogeneous  .i.. 
ture.   The  combustion  process  is  treated  as  a  constant 
pressure  combustion  and  immediately  followed  by  an  iscther- 

mal  mixing  between  newly  burned  an.^  ^k   u 

y    uurnea  and  the  burned  gases. 

As  proposed  by  Bracco  i-h^    ^ 

y  Bracco,  the  engine  combustion  models  can 

be  classified  as  follows: 

1.  Zero-dimensional  models 

2.  Quasi-dimensional  models 

3.  Multi-dimensional  models 

The  above  models  are  dxscu-ssed  earUer  in  Chapter  2 
2ero-dimensional  model  with  the  chosen  three  .ones 
control  volume  is  employed  n  this  analysis  because  it 
.ives  fairly  good  results  and  is  relatively  simple  to  use 
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4 . 2   Derivation  of  the  Energy  Related  burning  Velocity 

Equation 

The  burning  velocity  equation  is  derived  by  two  analy- 
tical methods:   Thermodynamics  analysis  and  combustion  me- 
chanism analysis. 


Assumptions : 

1.  During  combustion,  the  cylinder  contents  are  made 
up  of  three  distinctive  regions,  one  of  burned  part,  one  of 

unburned  part  and  the  flame  (reactive  region)  in  between 
separates  the  two  regions  above. 

2.  The  mass  burning  rate  of  the  fuel  air  mixture  is 
due  to  the  chemical  reaction  only  and  is  equal  to  the  rate 
of  mass  gained  to  the  control  volume  of  the  burned  part. 

3.  Pressure  surrounding  the  control  volume  of  burned 
region  is  uniform  throughout  the  combustion  process. 

4.  There  are  no  effect  of  the  potential  and  kinetic 
energy  in  the  combustion. 

5.  At  initial  stage,  there  is  a  very  small  (or  none) 
quatity  of  the  burned  mass. 

6.  The  burned  gases  consist  of  a  uniform  mi"ture  of 
products  of  combustion. 

7.  The  unburned  mixture  is  made  up  of  a  homogeneous 
mixture  of  air,  fuel  vapor  and  residual  products  of  com- 
bustion. 


JU 


8.  The  systems  behave  according  to  ideal  gas  law. 

9.  Heat  transfer  is  confined  to  convective  exchange 
between  the  gasses  and  the  walls. 

10.  The  energy  required  to  propagate  the  flame  from 
initial  to  final  position  is  denoted  as  the  flame  work  and 
is  defined  as: 


W     =W     -W      -W,     +W^ 
cv      exp    comp    loop    f 


(4.1) 


Where 


W   =  Total  net  work  for  a  given  system 
cv 

W     =  Expansion  work 
exp     ^ 

W     =  Compression  work 
comp 

W,     =  Loop  work 
loop       ^ 

W-;     =  Flame  work 


4.2.1   Thermodynamics  Analysis 

Fig. 4. 1  illustrates  the  initial  and  final  states  of  an 
arbitrary  chosen  control  volume  of  the  burned  gases.   Fig.  A-.Z 
shows  the  system  after  the  combustion  of  an  charge  element 
^m  and  the  uniform  pressure  surrounding  the  control  volume. 

The  conservation  of  mass^  together  with  First  law  of 
thermodynamics  is  applied  to  the  system  as  follows. 


*  Conservation  of  Mass 

At  any  instant  during  the  process,  the  continuity 
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equation  is  written  for  the  system  as: 


Pif'cv  ^^  *    i^"e^^  -   f^^i 


dt 


(4.2) 


For  a  period  of  time  t,  equation  (4.2)  becomes 


^^•2  -  '^l^cv  +  ^^e  -  ^"^i  =  ° 


(4.3) 


where 


•    • 


m.  ,  m-,  =   Mass  rate  of  the  burned  gases  at  initial 

stage  and  final  stage  respectively. 

rfi.,  r   =   Mass  burning  rate  enter  to  and  esit  from 
i'   e  ^ 

the  control  volume  respectively. 
From  the  assumptions  2  and  5, 


m^  =   0 

m^  =   0 

e 

after  arranging  and  let  m,  denotes  the  mass  burning  rate, 
equation  (4.3)  would  become: 


m-  =   m .  =  m, 
2      lb 


(4,4  ) 


*  First  Law  of  Thermodynamics 

Writting  the  First  law  for  the  uniform  state,  uniform 
flow  (USUF)  process,  at  any  instant  of  time  for  the  control 
volume  of  the  burned  part,  during  the  process,  we  have: 


Kv.'"-'  i-  '-v 
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Initial 


P. 

1 

T. 


■"'Ill 


Final 


Figure  4.1:   Initial  and  final  state  of  the  system. 


^i--^':<^ 
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Figure  4.2:   Uniform  pressure  surrounding  the  system 
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Q    +  Z  m. (h.  + 

^  c  V        11 


V 


m  (h   + 
e   e 


+  gz.  )  =   (  — )    +  w 

1        ,^  cv     cv 
at 

V  2 

—^-      +  gz/)  (4.5) 

7  -e 


Def f erentiate  equation  (4.1) 


w    =  w     -  w      -  W,     +  W- 
cv     exp    comp    loop    r 


4.6 


•vvJ^i:* 


jilt 


V'-lvN'-v 


Substitude  (4.6)  into  (4.5)  and  solve  for  W^ 


W^   =      Q^^   +  Zmih.    -    ( 


dE 
dt 


)  -    W  +    W  +    W, 

cv  exp  comp  loop 


(4.7) 


where 

W^  =   Rate  of  flame  work  done  by  the  system. 

6   =   Rate  of  heat  transfer  in  to  the  system, 
^cv 

Z  m.h.  =   Rate  of  total  enthalpy  input. 

(^)    =  Rate  of  energy  accumulated  in  the  system 
cv 
dt 

W      =   Rate  of  expansion  work  done  by  the  system, 
exp 

W      =   Rate  of  compression  work  done  on   the  system. 

comp 
• 
W,     =   Rate  of  loop  work. 

loop 

Let  dE-,/dt  denotes  the/^ight  hand  side  of  equation  (4.7) 
and  multiply  both  side  by  dt: 


cT  W^  =  m^^w^dt  =  dE^ 


(4.8) 


From  Fig.A.^,  pressure  surrounding  the  control  volume  is 
uniform,  the  work  done  is: 
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(p) (A) (dr^)  = 


^bO  "  ^bt 


=   F(dr. 


At  t  =  0;  r^^Q  =  0 


^f  =   ^^^bt^ 


Subsitude  (4.10)  into  (4.8) 


/w^  =  dE^  =  mj^(F)(rj^^)dt 


(F/A)  (A)  (dr^^) 


(4.9) 


(4.10) 


(4.11) 


From  definition  of  momentum,  per  Kg  mass  of  burned 
gases,  we  have: 


F  = 


dm, 

( )  (velocity ) /dm, 

dt 


velocity   _ 
dt 


u 


dt 


(4.12) 


Substitude  (4.12)  into  (4.11)  integrating  from  zero 
to  time  t  and  rearranging: 


Su  =   V^'^b^^^b^ 


^  =   W^/(m^)(r^) 


(4.13) 


Equation  (4.13)  represents  the  burning  velocity  as  a 
function  of  the  flame  work,  mass  burning  rate  and  the 
burned  part. 

This  equation  will  be  evaluated  term  by  term  and  the 


^y;<vv.' 
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mass    burning    rate: 


""fi  =  "'fo  ■"  "^""fl 


(4.14) 


Apply  the  First  law  to  the  flame  element , neglecting  the 
potential  and  kinetic  engergy,  we  have: 


dE 


cv    *•  1  1 


cv 


dt 


+  2ni  h   +  W 

e  e    cv 


and 


From  m  =  0 

e 

w=w    -w    -w,    +w^ 
cv     exp    comp    loop    f 


equation  (4.15),  after  solving  for  Vv'  ,  would  become 


(4.15) 


W^  =   Q    +  ^m^,h  ^  -   dE   /dt  -  W 
f     cv   ^    fl  uO      cv       exp 

+  W      +  W- 

comp     loop 


(4.16) 


Let  dE  /dt  denotes  the  right  hand  size  of  equation   4.16) 
and  from  definition  of  the  flame  work, 


"f  =  Pfi"^^fi'  =  ''fl 


and 


W 


w 


f  1 


f  1 


Pfi^^fi^^^i  -  ^0^  =  ^f^^i  -  ^0^ 


<f  w 


f  1 


(4.17) 


dt 


(4.18) 


Burning  velocity  is  defined  as  the  ratio  of: 
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computation  of  the  burning  velocity  will  be  given  in  sample 
calculation  in  Appendix  3. 

4.2.2   Combustion  Mechanism  Analysis 

Combustion  process  in  the  S.I.  engine  is  assumed  to 
progress  in  a  stepwise  maner  as  proposed  by  Patterson.   It 
is  described  below  with  the  assumption  of  no  mass  loss  due 
to  the  movement  of  exhaust  valve,  i.e  total  mass  contained 
in  the  cyclinder  remains  constant,  during  the  process. 

The  whole  system  consists  of  three  distinctive  regions, 
one  of  burned  part,  one  of  unburned  part  and  the  flame  be- 
tween separates  the  two  regions  above. 

The  combustion  is  considered  to  take  place  in  the  in- 
termediate processes  as  follows: 


Process 


0-1 


1-2 

2-3 
3-4 
4-5 


Constant  pressure  burning  of  an  increment  of 

the' charge  of  thickness  dr  =  r^  ~  ^o  ' 

Isothermal  mixing  between  newly  burned  and 

burned  gases. 

Isentropic  expansion  of  the  burned  gases. 

Constant  volume  heat  transfer. 

Polytropic  change  due  to  piston  motion. 


Table  I  summarizes  the  parameters  of  the  whole  system 
at  every  stage  of  the  process. 
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Initial 


nil, 


i!«l 


81 

j! 

-ml 


After  the  intake  and  isentropic  compression  processes 

in  the  cycle,  the  whole  system  is  assumed  to  have  a  negli- 
gible quatity  of  the  burned  mass.   Initial  thermodynamic 
properties  of  the  subsystems  are  given  as: 
Burned 


Mass 
Volume 
Temperature 
Pressure 


"bO  =  ° 

Not  applicable 
N.A 


Flame 


Mass 
Volume 
Temperature 
Pressure 


■"fO  =  ° 

Not  applicable 
N.A 


Unburned 


Mass 
Volume 


m 


uO 


uO 


Temperature 


uO 
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Pressure 


uO 


Step  1 ;   Constant  Pressure  Burning 

An  increment  of  mass,  thickness  of  dr  =  r.  -  r^  ,  is 
burned  at  constant  pressure.   The  flame  temperature  would  be 
at  its  adiabatic  flame  temperature.   Flame  element  will 
increase  its  volume,  thus  compressing  the  unburned  part. 
The  burned  gases  (if  any)  state  is  assumed  remain  unchange. 


The  whole  system  is  illustrated  in  the  above  figure. 
Let's  consider  each  of  the  subsystem  as  follows: 


Burned 


State  remained  unchange. 

Flame 

The  flame  properties  are  changing  from  the  initial 
condition  to  the  end  of  this  step  and  will  be  determined 
by  applying  the  First  law  to  the  flame  element.   Proper- 
ties at  this  stage  are: 
Pressure:    P^- 
Volume: 


'  .<.<■'. -Ij 
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S   = 

u 


dr 
dt 


^1  -  ^0 

dt 


(4.19) 


Substitude  (4.19)  into  (4.18)  and  solve  for  S^,  we 


have. 


S   =   W.,  /  F. 
u     r  I     r 


from   F^  =   m,- (a^  ^ ) 


(4.20) 
(4.21  ) 


and  acceleration  of  the  flame,  a^^,  is  defined  as: 


a    =   (Velocity  at  1  -  velocity  at  0)/dt 


^^1  -  ^0^ 
dt 


dt 


(4.22) 


At  initial  condition,  the  flame  radius  Tq  =  0. 
Substitude  (4.22)  into  (4.21)  then  (4.20)  and  noting  that 


^fl  ~  m-^  /  dt 


we  would  have: 


Su  =   W^  /  (m^^)(r^) 


(4.23) 


The  flame  work,  W^  will  be  determined  by  the  equation 
(4.16)  by  the  followed  processes. 


Step  2:    Isothermal  Mixing  Between  Newly  Burned  and 
Unburned  Gases. 
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A 

/ 

B 

/ 

u 

/ 

^1 

Burned  Gases 

The  volume  and  pressure  of  this  system  is  changing,  the 

final  volume  is: 

^b2  =   ^bl  ^^fl  =^fl                           "•'" 

Temperature  of  the  burned  gases  remained  constant  and 

is  equal  to  the  adiabatic  flame  temperature. 

'^b2  ^   "^fl  "  "^adia 

The  mass  of  the  system  is  equal  to  the  mass  entering            ' 

the  burned  gas.   Since  the  initial  mass  m^^  =  m^  ^  =  0                ,! 

was  assumed.   The  mass  of  the  system  after  the  mixing 

process  is. 

"^b2  =   "'fl  =  '"uO  -  ^u1                           ^""'^^^ 

Applying  the  equation  of  state  for  an  ideal  gas,  other 

properties  of  burned  system  can  be  calculated. 

Flame 

After 

the  isothermal  mixin 

g 

process,  the  flame  region 

'■■>-y.-'.'  ■ 
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is  reduced  in  its  size  with  a  negilible  thickness  and  is 
assumed  to  be  stationary  and  non  reactive,  until  the  end  of 
the  whole  process,  then  it  can  start  the  cycle  again. 


!i|!r 


Unburned  Mixture 

The  state  of  this  system  remained  unchange,  since  it  is 
assumed  that  the  mixing  process  occurred  between  the  burned 
and  newly  burned  systems  only. 

Step  3 :    Isentropic  Expansion  of  the  Burned  Gases. 


«■: 


Burned  Gases:   Expansion 

The  burned  system  undergoes  through  an  isentropic  ex- 
pansion process.   The  burned  mass  remained  constant.   With 
the  pressure  rapidly  increased  from  the  previous  process, 
the  burned  gases  would  compress  the  remalng  unburned  mix- 
ture and  causing  the  piston  to  move  downward. 

Properties  of  buned  gases  can  be  calculated  as  follow 


M. 


►      (K   -  1 )/K 
b3      p        P 


(4.26) 


b2 


b2 
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b3 


b2 


Where 


b2  ,   p 


K_  = 


b3 


:££. 


vp 


-  i  pi 

C   -  R 
PP 


(4.27) 


(4.28) 


The  expansion  work  can  be  determined  by  applying  the 
First  law: 


^exp  "   "b23 


=   Ub2^^b2)  -  "b3^\3) 


(4.29) 


Unburned  Mixture 

The  left-over  unburned  mixture,  in  constrast,  would  be 
compressed  by  the  burned  gas.   The  unburned  mass  is  reir.ained 
constant;  by  similiar  application  of  the  First  law  and  equa- 
tion of  state  for  an  ideal  mixture,  the  properties  of  this 
system  can  be  calculated. 

With  the  assumption  of  the  combustion  process  is  com- 
pleted, then  there  is  a  very  little  amount  of  unburned  T.ix- 
ture  presented. 


Step  4;    Constant  Volume  Heat  Transfer 

The  overall  heat  transfer  from  the  burned  to  unburned 
and  to  the  walls,  plus  heat  loss  from  unburned  to  the  wall 
can  be  computed  using  Woschni's  correlation.   Because  the 


'k   .  y 
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B 


'bu 


bw 


U 


"uw 


amount  of  heat  transferred  to  the  unburned  mixture  is  just 
off -set  by  that  from  the  unburned  gas  to  the  walls,  as 
proposed  by  Patterson.   The  heat  transferred  from  burned 
to  unburned  therefore  can  be  ignored,  the  convective  heat 
transfer  can  be  calculated  as  follow: 


and 


Qui  = 


h  (A  . ) (T  .  -  T  ) 
u   ui    ui     w 


u 


'bi 


=  ZQ 


(4.30) 


ui 


^b^^i^^^bi  -  ^w^ 


%i 


(4.31  ) 


Convective  heat  transfer  coefficients  are  computed 
from  Woschni's  correlation  and  is  given  as: 


.0.8 


h  =  no 


^0.2^^0. 53j 


Pi  1 


(4.32) 


The  meaning  of  each  term  in  this  equation  is  given  in 
Appendix  2. 
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The  wall  area  which  exposed  to  the  burned  gas  can  be 
calculated  as  given  in  the  below  equation: 


A  =   (TC)  Bore  [  stroke  ^^  ^  ^^^q. 


+  Rod( 1 


/   _   sin  Q   Stroke' 
y        4  Rod^ 


)  ] 


(4.33) 


The  Heat  transfer  area  for  burned  gases  is  calculated 
from  the  spherical  flame  front  model.   Area  for  the  unburned 
gases  is  the  difference  between  total  area  (of  chamber  walls, 
piston  and  head)  and  the  area  in  contact  with  the  burned 

gases . 

The  wall  temperature  is  measured  with  a  thermocouple, 
this  temperature  is  assumed  to  be  the  same  everywhere  on  the 
wall  surface  temperature  variations  are  quite  small  compared 
to  the  total  temperature  difference  between  the  burned  gases 
and  the  surfaces. 


m 


Step  5 :    Polytropic  Change  Due  to  Piston  Motion 


This  final  process  is  due  to  the  piston  motion,  where- 
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ther  it  is  polytropic  expansion  or  compression,  the  follow- 
ing relation  ship  can  be  applied  to  the  burned  and  unburned 
systems: 


V   k-1 


=   ( 


5  )  k 


(4.34) 


After  this  process  is  completed,  the  new  cycle  in  the 
cylinder  would  take  place  with  the  same  steps. 
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CHAPTER  V 


RESULTS  AND  DISCUSSION 


This  chapter  describes  the  results   obtained  from  the  theoretical 
and  experimental  work. 

Physical  and  Chemical  Properties 

When  evaluating  an   alternative  fuel,   knowledge  of   its  physical 

and  chemical  properties  is  essential  for  a  thorough  understanding  of 

how  the  proposed  fuel  should  perform  in  an  I.C.  engine.   While  there 

are  some   properties  which  are   generally  attributed  only   to  either 

S.I.   or  C.I.    engine  fuels,  many  properties  are  common  to  both,  as 

well  as  to  turbine  engine   fuels.   These  properties  include  density, 

specific  gravity,  absolute  viscosity,  surface  tension,   distillation 

temperatures,  flash  point,   and  heating  value.    The  fuel  properties 

which  have   the  greatest  effect  on  S.I.    engine  operation  are  the 

distillation   characteristics,   flash   point,    heating  value,    and 

composition.   In   the  present  work,    the  following  properties  were 

measured   or  computed:     1 .Di st i Hat  ion  characteristic,    2. Heating 

value,    3. Flash  point,    4. Water   tolerance,    5. API  gravity,    and 

6. Stoichiometric   air-to-fuel  ratio.    Most  of   the  properties   were 

measured  using  standard  ASTM  procedure  and  equipment. 

Distillation  Characteristics 

The  distillation  temperatures  for  the  15   fuels  in  this  study  are 
tabulated  in  Table  5.1.   The   data  for  indolene,   indolene-methanol , 
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and  indolene-MPHA  at  the  same  blend  levels  are  plotted  in 
Figs . 5.1-5. 5.  A  substantial  variation  in  distillation 
characteristics  is  evident.  Fig. 5.1  shows  the  distillation  curves  of 
pure  indolene,  methanol  and  MPHA.  Because  of  a  blend  of  many 
different  hydrocarbons  present  in  indolene,  it  shows  a  broad  range 
of  boiling  points.  Since  methanol  is  a  single  component  fuel,  it 
has  a  constant  boiling  point  temperature,  and  therefore,  does  not 
produce  a  distillation  curve. 

MPHA,  en  the  other  hand,  is  a  blend  of  different  alcohols,  and 
produces  a  curve  similar  to  that  of  indolene.  It  contains  75%  by 
volume  methanol,  and  a  remaining  25%  of  other  alcohols.  Therefore, 
the  distillation  curve  is  straight  line  for  about  75%  distilled  and 
then  varies  with  temperature.  The  shape  of  the  distillation  curve 
will  change  as  the  composition  of  MPHA.  For  example,  if  extra 
ethanol  or  other  higher  is  added  to  MPHA  before  mixing  with 
indolene,  the  distillation  curve  will  be  different  from  that  shown 
here . 

Figs. 5. 2-5. 5  show  the  distillation  characteristics  of  indolene, 
indolene-n-.rthanol  and  indolene-MPHA  fuels.  The  entire  distillation 
curve  is  affected  by  the  presence  of  alcohols.  The  constant 
temperature  boiling  due  to  methanol,  which  occurs  slightly  below  150 
T^  is  shown  clearly  in  all  figures.  At  all  blend  levels,  the 
indolene-MPHA  blend  is  closer  to  indolene  than  the  indolene-methanol 
blend. 

Although  distillation  is  a  continuous  process,   only  a  few  points 
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Fig.  5.1  Distillation  Properties  of  Pure  Fuels 
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Fig.  5.2  Distillation  Properties  of  20%  Blends, 
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Fig.  5.3  Distillation  Properties  of  40%  Blends. 
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Fig.  5.4  Distillation  Properties  of  60%  Blends. 
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Fig.  5.5  Distillation  Properties  of  80%  Blends. 
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along  the  curve  are  necessary  to  characterize  a  fuel,  namely,  the 
initial  boiling  point,  10  percent  maximum,  50  percent  ir.inimum,  50 
percent  maximum,  90  percent  maximum,  and  the  end  point.  Depending 
upon  the  volatility  class,  the  maximum  10  percent  limit  varies  from 
122F  to  158F,  minimum  50  percent  is  170F,  maximum  50  percent  varies 
from  230F  to  250F,  maximum  90  percent  varies  from  365F  to  374F  and 
end  point  is  437F.  Since  the  50  percent  minimum  limit  is  20  to  25 
degrees  above  the  constant  boiling  point  temperature,  that  part  of 
the  distillation  specification  will  be  difficult  to  meet  using 
indolene-methanol  blends.  As  shown  in  Fig. 5. 2,  the  20  percent 
alcohol  blend  is  within  the  specified  limit.  The  40  percent  blend, 
Fig. 5. 3,  exceeds  the  limit.  However,  by  modifying  the  formulation 
of  MPHA,  it  may  be  be  possible  to  meet  this  part  of  the 
specification  upto  70  percent  blend  level. 


The  shape  and  level  of  the  distillation  curve  affects  gasoline 
engine  performance.  The  characteristic  indicating  the  cold  starting 
quality  of  a  fuel  is  represented  by  0  to  20  percent  evapordted.  A 
highly  volatile  0  to  20  percent  fuel  will  aid  cold  starting,  but  may 
tend  to  cause  vapor  lock  at  high  teperatures.  The  warm-up 
characteristics  are  related  to  the  temperature  at  which  50  percent 
of  the  fuel  is  vaporized.  A  low  50  percent  point  yields  better 
performance,  but  with  an  increase  in  the  tendency  to  cause 
carburetor  icing.  Fuel  economy  relates  to  heating  value,  and 
density,  which  may  also  correlate  somewhat  with  the  90  percent 
point.  Low  back-end  volatility  may  lead  to  crankcase  dilution, 
deposits,  gum,  and  spark  plug  fouling. 


56 


H^H 

"ll 

'-^H 

lit 

"^1 

.1" 

.^l^^l 

II 

''^^911 

Ik 

:H 

Ml 

^■■•^ 

•'•'  ■'■'■""•' 

,,fi 

"c;^;--.-^:-- 

^>;i;;>A' 

-I'll! 

^'■'  .^s' 

ii. 

'Ill 

iirj 

Energy  Density 

Fuel  economy  is  usually  expressed  in  terms  of  BSFC  for  engines. 
The  BSFC  is  a  function  of  the  energy-mass  density,  which  affects  the 
conversion  effiency  of  the  chemical  energy  to  useful  mechanical 
work.  The  energy  density,  or  heating  value  of  a  fuel  is  a  measure  of 
the  amount  of  energy  released  as  a  result  of  the  complete  combustion 
of  the  fuel.  The  higher  heating  value  (water  formed  by  combustion  is 
in  vapor  form)  is  measured  using  an  adiabatic  bomb  calorimeter.  The 
heating  values  of  indolene,  indolene-methanol ,  and  indolene-MPHA 
blends  were  measured.  Fig. 5.6  shows  the  heating  values  of  the 
indolene-alcohol  blends.  The  measured  heating  value  decreases  as  the 
concentration  of  alcohol  increases  in  the  blend.  The  heating  value 
of  indolene-MPHA  blend  is  higher  than  indolene-methanol  blend  at  all 
blend  levels. 

The  energy-volume  density,  BTU/gal,  is  dirctly  proportional  to 
the  produst  of  the  energy-mass  density  and  the  specific  gravity. 
Since  the  specific  gravity  of  MPHA  is  similar  to  that  of  methanol, 
the  energy-volume  density  of  indolene-MPHA  blends  also  be  great-er 
than  that  of  indolene-methanol  blends. 


The  heating  values  of  blends  are  generally  predicted  from  the 
pure  component  properties  by  equation  5.1.  Figure  5.7  shc'S  the 
blend  heating  values  predicted  from  its  pure  components  using  this 
equation,  end  the  measured  values,  for  indclene-MPHA  blends. 

KVM  ^  MI  *  HVI  -r  MA  *  HVA  (5.1) 

where  HVM  =  Heating  value  of  mixture 
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Fig.  5.6  Measured  Higher  Heating  Values  of  Indolene-Alcohols  Bleni 
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MEASURED  i  CALCULATED  HEATING  VALUE 
INDOLENE-MPHA  BLENDS 
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MI  =  Mass  fraction  of  indolene 
HVI  =  Heating  value  of  indolene 
MA  =  Mass  fraction  of  alcohol 
HVA  =  Heating  value  of  alcohol 

There  is  a   good  agreement  between  measured  and  predicted  values 
for  indolene-MPHA  blends. 


Water  Tolerance 

Fig. 5. 8  shows  the  water  tolerance  of  the  indolene-alcohol  blends. 
The  solubility  of  alcohols  in  indolene  when  water  is  present  even  in 
minute  quantities  is  limited  at  room  temperature.  When  small  amounts 
of  water  added  to  blends  of  alcohol  with  indolene,  hydrogen  bonds 
form  between  the  water  and  alcohol  molecules,  and  the  blends 
separate  into  two  phases.  Paraffinic  hydrocarbons  predominate  in  the 
upper  phase,  while  the  lower  phase  consists  primarily  of  alcohols 
and  water.  Water  tolerance  is  also  significantly  decreased  at  lower 
temperatures.  At  room  temperature,  as  the  alcohol  concentration 
increases,  water  tolerance  increases.  The  90  percent  methanol  blend 
has  a  water  tolerance  of  8.25  percent  by  volume,  compared  to  0.25 
percent  at  10  percent  blend.  The  90  percent  MPHA  blend  has  about  50 
percent  higher  water  tolerance  compared  to  the  90  percent  methanol 
blend  at  room  temperature. 

Solubility  of  methanol  and  indolene  in  the  presence  of  water  can 
be  one  of  the  most  difficult  problems  attending  the  use  of  these 
blends  as  motor  fuels.  The  fuel  distribution  system,  refinery, 
pipeline,   and  terminal  tanks  now   used  are  not  necessarily  designed 
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Fig.  5.8  Water  Tolerance 
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to  fully  exclude   rainwater  from  coming  into  contact   with  the  fuel. 

Condensation   in   the   storage   tanks  can   cause   additional   water 

contamination.   Using   MPHA-gasoline  blends   instead   of   methanol- 

gasoline   blends  increases   the  water  tolerance,   and  may  be   used 

without  making  major  modifications  to  the  existing  fuel  distribution 
and  delivery  system. 

Flash  Point 


The  flash  point  is  the  lowest  temperature  at  which  the  vapors 
arising  from  a  fuel  will  momentarily  ignite  as  a  flame  is  swept 
across  its  surface.  It  is  used  to  assess  the  flammability  of  a  fuel 
and  to  establish  shipping  and  safety  regulations.  The  flash  points 
of  methanol  and  MPHA  fuels  were  measured  using  a  Pensky-Martens 
closed  up  tester,  following  the  ASTM  D93  method.  The  flash  point  of 
pure  methanol  was  50F,  and  that  of  MPHA,  49F. 

API  Gravity  and  Specific  Gravity 

The  API  gravity  is  defined  by  an  arbitrary  scale  whose  units  are 
degrees  API,  and  is  measured  at  15.5  deg.C.  The  API  gravity  of 
indolene-methanol  and  indclene-MPHA  blends  were  measured  using 
Fisher  brand  Petroleum  hydrometers.  Specifications  for  API  gravity 
hydrometers  have  been  established  by  ASTM.  Fig. 5.9  shows  the  API 
gravity  of  indolene-methanol  and  indolene-MPHA  blends.  In  genera], 
as  the  alcohol  content  increases,  API  gravity  decreases.  Indolene- 
methanol  has  a  slighcly  higher  API  gravity  than  indolene-MPHA. 

The  specific  gravity  is  found  by  dividing  the  density  of  the  fuel 
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Fig.    5.9    API    Gravity 
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by  the   density  of  water.    The  specific   gravity  is  related   to  API 
gravity  as  shown  in  equation  5.2. 

Degree  API  =  ( 141 . 5/spec i f ic  gravity)  -  131.5        (5.2) 

Fig. 5. 10  shows  the  specific  gravity  of  indolene-methanol  and 
indolene-MPHA  blends.  As  the  alcohol  concentration  increases,  the 
specific  gravity  increases  slightly.  In  general,  indolene-MPHA  has  a 
higher  spefic  gravity  than  indolene-methanol  blend. 

Stoichiometry 

Stoichiometry  and  combustion  are  related  to  the  H/C  ratio  of  the 
fuel.  In  general,  for  indolene,  H/C  is  1.86,  4.0  for  methanol,  and 
3.269  for  MPHA.  The  stoichiometric  air-fuel  ratio  depends  upon  the 
H/C  ratio  of  the  fuel.  The  stoichiometric  air-fuel  ratio  is  14.56 
for  pure  indolene,  6.43  for  methanol,  and  7.4468  for  MPHA.  Fig. 5.11 
shows  the  variation  of  stoichiometric  air-fuel  ratio  as  a  function 
of  alcohol  concentration.  As  the  alcohol  concentration  Increases, 
stoichiometric  air-fuel  ratio  decreases.  The  stoichiometric  air-fuel 
ratio  of  indolene-MPHA  blends  is  higher  than  those  of  indolene- 
methanol  blends. 


PERFORMANCE  PARAMETERS 

Effect  of  Alcohol  on  MBT  Spark  Advance 

Figures  5.13-5.16  show  the  effect  of  varying  the  methanol  and 
MPHA  concentrations  on  the  MBT  spark  advance,  at  different 
compression  ratios.  At  any  fixed  compression  ratio,  as  the  alcohol 
concentration  increases,  the  MBT  spark  advance  decreases,   or  shifts 
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Fig.  5.10  Specific  Gravity 
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MBT  SPARK  ADVANCE  AT  CR=5.0 
INDO.-METH.  &  INDO.-MPHA  BLENDS 


I'll 


III 


o 

Q 

H 

m 
d) 

■D 

UJ 

O 


D 
< 

DC 

< 

Ql 


/  u.u  - 

1 

i 
1 

1 
1 

LEGEND 

1 

1       i 

■  1 
j             i 

ftn  n 

•  =INDO.+METH. 

i 

! 

i             ' 

:        i 

DU.U 

L 

■  =1 

MDO. 

+MP 

HA 

1 

i 

j 

i 

1 

i 

i 
! 

50.0- 

■^ 

i               ~ 

1^^^ 

..  .  1. .       

! 

1 
1 
! 

1 

1 
1 
1 

AH  n 

^ 

1 

•rU.U 

w   ^^^^^ 

30.0- 

1 

. 

! 

20.0- 

1 

i 

1 

0.0     0.1 


0.2      0.3      0.4      0.5      0.6      0.7      0.8      0.9      1.0 

FRACTION  OF  ALCOHOL 


Fig.    5.13   MBT  Spark   Advance   at   CR=5:1 
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Fig.    5.14   MBT   Spark   Advance    at   CR=6:1 
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Fig.    5.15   MBT   Spark  Advance   at   CR=7:1 
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Fig.  5.16  MBT  Spark  Advance  at  CR=KLCR 
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towards  the  top  dead  center.  This  is  an  indication  of  a  higher 
burning  velocity  of  indolene-alcohol  blends  compared  to  purej 
indolene.  At  any  compression  ratio,  the  MBT  spark  advance  of 
indolene-MPHA  blends  is  higher  than  that  of  indolene-methanol 
blends.  The  MBT  spark  advance  is  controlled  by  the  thermodynamic 
properties,  such  as  temperature,  pressure  and  burning  velocity  at 
one  specific  engine  operating  condition.  With  the  same  intake  air 
temperature,  the  temperature  of  the  air-fuel  mixture  during  the 
compression  process  is  related  to  the  latent  heat  of  the  fuel.  MPHA 
has  a  lower  latent  heat,  and  a  higher  burning  velocity  than; 
methanol.  However,  the  combined  effect  on  spark  advance  is  a  higher 
MBT  spark  advance  in  the  case  of  indolene-MPHA  blends. 

Effect  of  Alcohol  on  KLCR 


Table  5.2  shows  the  KLCR  as  a  function  of  methanol  and  MPHA 
concentration  in  the  blend.  As  the  methanol  concentration  increases 
from  0  to  40  percent,  KLCR  increases  from  7:1  to  10:1.  Further 
increase  in  methanol  concentration  does  not  cause  any  further 
increase  in  the  KLCR.  The  variation  of  KLCR  with  MPHA  fraction  is 
similar  to  that  of  methanol.  However,  the  knocking  intensity  with 
indolene-MPHA  was  lower  than  indolene-methanol  blends.  The  maximum 
compression  ratio  that  can  be  obtained  in  this  type  of  engine  is 
10:1.  The  KLCR  of  indolene-MPHA  is  higher  than  indolene-methanol 
blend.  The  actual  magnitude  of  the  KLCR  could  not  be  obtained, 
however,  due  to  the  limitation  of  the  engine. 
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TABLE  5.2 
KNOCK  LIMITED  COMPRESSION  RATIO 


Volume  Fx act ion  of 

Knock  Limited 

Alcohol  in  Indolene 

Compression  Ratio 

0.3  aethanol 

8.5 

0.5  Methanol 

10 

0.7  aethanol 

10 

1.0  nethanol 

10 

0.3  MPHA 

8.87 

0.5  MPUA 

10 

0.7  MPUA 

- 

1.0  MPUA 

^ 
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Effect  of  Alcohol  on  Brake  Power 

Figures  5.17-5.20  show  the  brake  power   as  a  function  of  fraction 
of  alcohols.   In  general,   as  the  alcohol   concentration  increases, 
brake  power  increases.   At  any  compression  ratio,  or  any  fraction  of 
alcohol,   the  brake  power  of   indolene-methanol  is  no  different  from 
that  of  indolene-MPHA. 

Effect  of  Alcohol  on  Brake  Specific  Fuel  Consumption 

Figures  5.21-5.24  show  the  brake  specific  fuel  consumption  as  a 
function  of  fraction  of  alcohols.  In  general,  as  the  alcohol 
concentration  increases,  the  brake  specific  fuel  consumption 
increases.  At  any  compression  ratio,  indolene-methanol  has  a  higher 
BSFC  compared  to  indolene-MPHA,  specifically  after  70  percent  blend 
level.  It  is  noted  that  BSFC  is  on  a  mass  basis,  (kg/kW-hr).  With 
regards  to  BSFC,  methanol  is  at  a  distinct  disadvantage  compared  to 
MPHA  because  of  its  lower  heating  value. 

Effect  of  Alcohol  on  Brake  Thermal  Efficiency 

Figures  5.25-5.28  show  the  brake  thermal  efficiency  as  a  function 
of  alcohol  fraction.  At  any  compression  ratio,  the  thermal 
efficiency  of  indolene-alcohol  blend  is  higher  than  that  of 
indolene.  In  general,  the  indolene-methanol  blends  have  a  higher 
thermal  efficiency  than  indolene-MPHA  blends.  The  brake  thermal 
efficiency  depends  upon  combustion  period,  cooling  losses,  chemical 
equilibrium  losses,  change  in  specific  heat  losses,  and  latent  heat 
of  vaporization.   The  thermodynamic  analysis  of  the  entire  cycle  is 
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Fig.    5.17   Brake   Power   at  CR=5:1 
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Fig.    5.18   Brake   Power   at   CR=6 : 1 
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Fig.    5.19    Brake   Power    at   CR=7:1 
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Fig.    5.20   Brake   Power  at   CR=KLCR 
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Fig.    5.21   BSFC   at  CR=5:1 
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Fig.     5.22    BSFC    at   CR=    6:1 


^^^i^ 


79 


O 

LL 
CO 

CD 


1.2 


1.0 


0.0 


FUEL  CONSUMPTION  AT  C.R  =7.0 
INDO.-METH.  &  INDO.-MPHA  BLENDS 


LEGEND 

•  =  INDO.-METH. 
■  =  INDO.-MPHA 


;>:V-;'k;Jv 


0.2  0.4  0.6  0.8 

FRACTION  OF  ALCOHOL 


Fig.    5.23    BSFC   at   CR=7:1 


80 


1.2 


1.0 


O 
CO 


FUEL  CONSUMPTION  AT  KLCR 
INDO.-METH.  &  INDO.-MPHA  BLENDS 


0.0 


LEGEND 

=  INDO.-METH. 
=  INDO.-MPHA 


0.2  0.4  0.6  0.8 

FRACTION  OF  ALCOHOL 
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Fig.  5.24  Thermal  Efficiency  at  CR=5:1 
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Fig.    5.26    Thermal    Efficiency   at   CR=6:1 
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Fig.  5.27  thermal  Efficiency  at  CR=7:1 
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Fig.  5.28  Thermal  Efficiency  at  CR=KLCR 
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■    required  to  explain  the  effect  of  indolene-MPHA  on  thermal  effiency. 

Effect  of  Alcohol  on  Cycle-to  Cycle  Variation 

Fig. 5. 29  shows  the  maximum  standard  deviation  of  the  cylinder 
pressure  in  a  cycle,  as  a  function  of  fraction  of  alcohol,  and  also 
as  a  function  of  compression  ratio.  In  general,  as  the  alcohol 
concentration  increases,  the  standard  deviation  increases.  In  other 
words,  alcohol  increases  the  cycle-to-cycle  variation,  although  only 
marginally.  The  factors  that  affect  cycle  to  cycle  variation  could 
be  spark  scatter,  inhomogenei t ies  in  mixture  temperature,  ignition 
delay,  and/or  fluid  motion  in  the  vicinity  of  the  spark  plug. 
Howwever,  a  more  detailed  analysis  of  the  ignition  delay, 
turbulence,  and  other  parameters  would  have  to  be  made  before  the 
exact  causes  for  the  increased  cycle-to-cycle  variation  can  be 
ascertained. 

Effect  of  Alcohol  on  Flame  Speed 

Figures  5.30-5.40  show  experimentally  measured  flame  speed  as  a 
function  of  crank  angle  during  the  combustion  process.  No  clear 
trend  can  be  observed  as  to  the  effect  of  alcohol.  This  can  be 
explained  as  follows: 


The  flame  speed  is  function  of  the  pressure  and  temperature  of 
the  unburned  gases,  turbulence,  equivalence  ratio,  and  the  fuel 
type.  Since  the  intake  air  was  kept  constant  during  this  study,  the 
mixture  temperature  varied  according  to  fuel  blend  used.  As  the 
alcohol  content  increased,   mixture  temperature  decreased,   lowering 
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the  flame  speed.  Had  the  experiments  been  conducted  at  a  constant 
mixture  temperature,  a  higher  flame  speed  would  have  been  obtained 
with  the  alcohol  blends. 


Effect  of  Alcohol  on  Burning  Velocity 

As  with  the  flame  speed,  the  burning  velocity  is  also  affected  by 
the  mixture  temperature.  It  is  therefore  inappropriate  to  compare 
the  burning  velocities  of  indolene-alcohol  blends  with  those  of 
indolene.  It  can  be  seen  in  Figures  5.41-5.44  that  at  any  fraction 
of  alcohol,  indolene-MPHA  blends  have  a  higher  burning  velocity 
compared  to  indolene-methanol  blends.  However,  the  difference  is 
only  marginal. 
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Fig.  5.41  Burning  Velocity  of  10%  Blends. 
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Fig.  5.42  Burning  Velocity  of  20%  Blends 
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Fig.  5.43  Burning  Velocity  of  30%  Blends 
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Fig.  5.4  4  Burning  Velocity  of  40%  Blends. 
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CHAPTER  VI 

CONCLUSIONS 

The   following   are   the   conclusions   of   the   theoretical   and 
experimental  work. 

1.  The  blend  containing  more  than  20  percent  alcohol  did  not  meet  50 
percent  minimum  distillation  specification. 

2.  The  heating  value  as  measured  with  adiabatic  bomb  calorimeter,  of 
indolene-MPHA  is  higher  than  that  of  indolene-methanol  blend  at 
all  blend  levels. 

3.  The  water  tolerance  of  indolene-MPHA  blend  is  higher  than 
indolene-methanol  blend  at  all  blend  levels. 


4.  The  API  gravity   of  indolene-methanol  is  higher   and  the  specific 
gravity  is  lower  than  indolene-MPHA  blend. 

5.  The  stoichiometric  air-fuel  ratio  of  indolene-MPHA  is  higher  than 
indolene-methanol  blend. 

6.  Adding  methanol  or  MPHA  to  indolene   causes  the  KLCR  to  increase. 
The  KLCR  of  MPHA  is  greater  than  that  of  methanol. 

7.  MBT   spark   timing   increases   as   the   alcohol   concentration 
.^    increases.   Indolene-methanol  has  a  lower  MBT  spark  advance  than 

indolene-MPHA. 

8.  As  methanol  and  MPHA  concentration  increases  in  the  blend,   power 
output,  BSFC,   and  brake  thermal  efficiency  increase.    Indolene- 
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MPHA   has   the   same   power  output,    lower   BSFC,    and   thermal 
efficiency  compared  to  indolene-methanol  blend. 

9.  The  actual  burning  velocity  of  indolene-MPHA  are  slightly 
higher  than  indolene-methanol  blends  for  the  blends  upto 
40  percents.  However  there  are  mixed  trends  for  flame  speed. 
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